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^ _ Abstract 

> [ 

'nI" ! We consider the generalized Ornstein- Uhlenbeck equation dtX = —mXt + rj. In this paper We 
construct the Levy noise rj. The generalized Ornstein- Uhlenbeck process Xt will be represented 

■ by a special types of graphs called rooted trees with two types of leaves. 
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X : 1 Introduction 



In this paper we study the generalized Ornstein- Uhlenbeck (OU) equation: 

§lXt{x) = -mXt{x) + ri{t) , m > . . 

Xo{x) = 0(x), {t,^{x)) G]0, oo[yi^'^ ^ ' 

Where rj = r]{t) is a Levy noise (more details will be given in section 2). 

Generalized Ornstein-Uhlenbeck processes has been introduced by Barndorff-Nielsen (1998). 
However by replacing the continuous time i by a discrete one, we introduce the generalized 
random process, then by the use of the well known Bochner-Minlos theorem, see eg. 0, we 
construct the Levy noise r]. 

A graphical calculus will be applied to the generalized OU-process in this paper. In fact we will 
introduce a special types of graphs called rooted trees with two types of leaves, see. eg [U |5] an 
analytic value will be given to each rooted trees and therefore a graphical representation of the 
OU-process will be recalled in this paper. 
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The procedure of given each graph a numerical value is called "perturbation theory". However 
our model of representing the solution of equation ([T|) in terms of rooted trees with two types of 
leaves is not restricted to such equation in fact one can generalize equation ([1]) and make it more 
complicated, e.g by introduction of non linear terms, e.g for force F of gradient type, F = W 
we obtain a non-linear SDE and one can ask the same questions as before, again expansion into 
graphs is possible. Generalizing further we pass from SDE's driven by Levy noise to SPDE's. 
Graphical representation of the OU process seems to be of a great importance in fact they can be 
used for modeling biological processes such as neuronal response, see e.g [7j, also in Mathematical 
Finance, see. eg [6j, the modeling of the dynamics of interest rates and volatilities of asset prices 
where each notions will have a specific graph representation. 
The remainder of the article is organized as follows: 

The next section will be reserved to the construction of Levy noise, we will develop a interesting 

model which help the reader to understand the other sections of our paper. 

In section 3 we recall a theorem which determine the solution of equation ([1]), the distribution 

of the generalized OU process will be given as well as a particular case when rj will be Gaussian 

noise. 

In section 4 we introduce the rooted trees with two types of leaves. We will develop an algorithm 
which gives a numerical value to each rooted tree, in the last part of this section we will recall 
a graphical representation of the generalized OU process. 

2 Construction of the Levy noise 

This section is devoted to the construction of the Levy noise rj, so far we introduce the Levy 
characteristic ip. By taking the time to be discrete the Fourier transform of the joint distribution 
of r] will converge to a function which depend of ip, the generalized random process will be given 
as well. Hence by the use of the well known Bochner Minlos theorem we recall a theorem which 
gives the Levy noise. 

It is well known that an infinite divisible probability distribution P is a probability distribution 
which satisfy the existence of a probability distribution such that P = P^* - ■ ■*Pn (n times). 
Let P be an infinitely divisible probability distribution. By Levy-Khinchine theorem (see e.g. 
[3]), the Fourier transform, or characteristic function, of P satisfies: 

1. Cp{t) = /jg e*^'^'*^ (iP(,s) = e^^^\ t £ R, where ■0 is the Levy characteristic of P which is 
uniquely represented by ,see. eg Lukacs [12], 

i){t) = iat + z (e*"* - 1) dM{s), V t G M. (2) 

2 Jr\{o} 

with M satisfies /ig\^|Q} min(l, s'^)dM{s) < oo. 
We can explain the meaning of ([2]) as follows: 

If 0-2 = M = 0, then ipit) = iat and the process {X{t), t > 0} is simply deterministic motion, 
here a is called the drift and it's physical interpretation is the velocity of this motion. 
Now if only M = so that = iat - \aH'^ and therefore Cp{t) = el'"*"!'^^*^] this is a 
characteristic function of a Gaussian random variable X{t) having mean ta and covariance ta. 
The process t > 0} in this case is known as the Brownian motion. 

The last case is when a 7^ 0, o" 7^ and M 7^ 0, in this case the process {X{t), t > 0} can be 
represented by X{t) = bt + ^/aB{t) + N(t) where h is the drift, B{t) the Brownian motion and 
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N{t) will be recognized as a poisson process with intensity A taking values in {nh n G N}. 

In the following we denote by F the probability distribution on of a given stochastic process. 

Definition 2.1. A Generalized stochastic process r](t) is a stochastic process such that the fol- 
lowing proprieties are satisfied: 

. F{r]{t)) = F{r]is)) , \/t,s>0 

• ri{t) is independent of r]{s), i.e F{r]{t)r]{s)) = F{r]{t)) (g) F(r/(s)) 

From now our stochastic process rj will be a generalized stochastic process and we would like 
to construct a generalized Levy noise: 

Let us start with a model, where the continuous time t G M+ is replaced by discrete time 
t G -N. In this case one can thus model the noise {'n^^\t)}t^ J-n ^ collection of i.i.d. random 

n 

variables such that the "global random fluctuations" that the noise gives should not depend on 
the lattice scale ^ fixed in our model, i.e to model the noise independently of the lattice scale, 
we need that: 

F(r?«(l)) = Ffr?(")(-) + 7?(")(-) + ---+r?(")(-)' 
\ n n n ^ 

V(r?W(l)))*" (3) 
Here F{ri'^^\l)) is the random fluctuation of the unit interval. 

Definition 2.2. The Fourier transform of rj^^^t) is given by 

-^(r/("^)(/)(i) = e^,/(t) G M'^ (4) 

Where ip is the Levy characteristic given in equation 

Proposition 2.3. The Fourier transform of the joint distribution of r]^"\^), ??''"'^(^), •.• con- 
fa ■>P{f{i)) dt 

verges to e + , as n — > oo. 



Proof, we have: 



T if i-).f {-)..:) =11^"^ 

\ n n J n n 



E 



e 



i'(f(t)) 



T ■ 111 ■ /e ^{fW)<it 

It IS now easy to see that the last expression converges to e + , as n — > oo. 

■ 

Now it turns out that, though the above limit exists, it does not give us the definition of rj(t) as 
a random function of as the limit: 



F 



((r?W(l))) =5o (5) 



shows that the fluctuations of the random variable ri^"'\t) vanish for n — > oo. 

The solution of this apparent paradox lies in the fact, that r]{t) is not a stochastic process in the 

classical sense, but rather a distribution in t then a function of t. 
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Definition 2.4. A generalized random process or a random field is a mapping: 

4> : {n, B, P) 5'(M, M'^) 

such thatoj — > {(t){w), f) is measurable and for fn — > f inS{R, W^), F{(f){fn)) — > F{(f){f)) 
in low, i.e. 

lim / F{<p{U))dP = lim / F{x)F{<p{Umx)dx 
= f F{x)F{ct>{f)){x)dx 

= [ F{m)dP 
Jn 

Theorem 2.5. (P. Levy) Let X and X„, n e N be a stochastic processes then: 
Xn ^ X in low if and only if F{Xn){k) — > J^{X){k), eW^. 

Definition 2.6. Two random fields (pi and 02 are equivalent in low if all finite dimensional 

distributions of (f)i and (f)2 coincide, 

i.e. //V/i,...,/n G M'^) and'iAu.-.^An e H(M) 

P{Mh)^-Al,..;Mfn) G X) =P(02(/l) e A,...,02(/n) £ An) 

which means that: 

F{Mh) ■ ■ ■ Mfn)) = F{Mfi) ■ ■ ■ Mfn)) (6) 

But how can one construct such random field? 
We note by B the cr-algebra generated by the cyhnder sets of 5'(M'^). Then (5'(]R'^), B) is a 
measurable space. 

We define a characteristic functional on iS(M'^), as a functional C : iS(M'') — > C such that : 

1. C{0) = 1. 

2. C is continuous on S'(M'^); 

3. C is positive-definite, i.e. Vzi, Zn G C,Vn G N, /i, /„ G S{R^) 

n 

^ zizjC{fi-fj) >0. (7) 
i,j=o 

In the following we choose (Q, B) = {S' , B{S')). 

By the well known Bochner Minlos theorem, we know that for a given a Characteristic functional 
C, there exist a unique (up to equivalence in low) random field (f) such that: 



C{f) = Efe^-^^^)] = / e^'^(-^)H dP{uj), V / e S{R, R'^) 
Jn 



(8) 
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Theorem 2.7. Let ip be a Levy characteristic given by the representation ^ then there exist a 
unique measure rj on S' such that: 

C{f)= [ e*<'^'-^>d?7(w) =exp ( / tl;{f{x))dx), x G M^. (9) 

Js' yJRd. J 

Proof. It is easy to see that the right hand side of Q is a characteristic equation, now Bochner 
Minlos theorem concludes. 



Definition 2.8. The measure r] given by theorem \2. 7| j is called white measure with Levy char- 
acteristic il), and {S'(R'^), (3, rf) will be the generalized white noise space associated with ip. The 
associated coordinate process 

7? : S{R'') X (S'(r), B, ri) C, = (/, V/ G 5(R'^), ^ G S\T) (10) 

is called a Levy noise. 



3 Generalized Ornstein-Uhlenbeck processes 

In this section the generahzed stochastic process rj wih be taken as a Levy noise in the sense 
of the previous section. We recall a result which gives the distribution of the generalized OU 
process, this later is known as "Mehler's Formula." 
When there is no confusion we note X{t) by Xf. 

Definition 3.1. The stochastic process {Xt)t>Q verifying equation (OP, where rj is a Levy noise 
is called the generalized Ornstein- Uhlenbek (OU) process. 

Theorem 3.2. The solution of the generalized OU- equation ([ip is given by: 

Xt{x) = e-"^' <i){x) + Te"'" r/(t')di' (11) 







Proof. The solution of the Homogeneous equation associated to equation ([TJ is given by 
Xt^c = C e~™*, now by the method of variation of parameters one thus get 



t 



Xt{x)= / e-™^*-*^7?(t')dt'+ Ce-"* (12) 
Jo 

The result now follows by taking Xq[x) = C = (j){x). 



The problem now is how to define the integral in theorem ()3.2p ? 
Let Cj = (0, 1, 0...) G W^, we put 

xL,*(0 = e,-l*'<ie— (13) 

and 

X^t) = e-^' cP^ + v{xL,t) (14) 

This looks better, but Xm t ^ ^^^^ function in ^(IR, W^). 

Let 

Xk,„(t') = e,x"(i-0e-™^*-*'^ (15) 
Here x" is a sequence of C°°— test functions that approximate the function from below. 
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Proposition 3.3. The sequences {Xn}n>i converges in P), as n — > oo. 

Proof. 



E 



E 



V [Xm,t,n Xm,t,n' 



where a, /3 are constants, now by the use of the dominated convergence theorem it is clear that 
the last expression goes to as n, n goes to infinity. 



Let F be the distribution of the generalized OU process, then the following result holds: 
Theorem 3.4. The distribution F of the generalized OU process is given by: 

F(Xt) = F(e— %)*F(r?(x„.,t)) 
Proof. The Fourier transform of Xt is given by: 



(16) 



= exp (^i {Xo, p)e-'^' + ^{e~^^'-'">p)dt'^ 



Corollary 3.5. 1)- Let r] be a Gaussian noise, then the distribution of the OU process is given 
by: 



F(Xt)(x) = (detA)-i(— (l-e-^-*))-! exp.^ ^ 



m 



2 1 - 



where A = (£(1 - e"^™^)) \ 
2)- The distribution of the Brownian motion is given by: 



(17) 



F{Bt){x) = (cZeU)-^(47r)-i exp |-l((x - Xq), D-\x-Xo))jdx 

Proof. 1)- If the noise rj is taken to be Gaussian noise, then ipi^) = —{k, Dk), thus: 

r V;(e-'"(*-*')/c)di' = -{k, Dk) f e-^'^^^-^'hdt' 
Jq Jo 



(18) 
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Hence by applying the Fourier integral, one thus obtains: 



VdetA f 1 
2 



F(Xi)(x) = exp --{{x - e-"*Xo), Aix - e-™*Xo)) }dx 



(2vr)- 

= (deM)-i(^(l-e-2-*))-f exp|-i^— ^^((x-e— *Xo),D-nx-e— *Xo))}dx 
Here^= (£(l-e-2™*))-\ 

2)- If m — > the right hand side of equation (I17p converges to 

(cletAy^ATr)-^ exp - Xq), D-\x - Xo))^dx (19) 

which the distribution of the Brownian motion. 



Remark 3.6. 1. The expression ofF{Xt) given by equation (11) is known as "Mehler's For- 
mula". 

2. For the Gaussian case, i.e, ip{k) = —{k,Dk) one thus obtain the Heat equation: 



Po(x) = f{x) 

3. The formula given by equation < f7?|) is the semigroup of the Brownian motion. Such a formula 
for the non Gaussian case is missing ! 

If now 'ilj{k) is taken to be of the more general type, i.e, ip{k) is a Levy-Khinchine function, 
see e.g [3], one thus get: 



J=l J j,l=l J ' 



+ z [Pt{x + y)-Pt{x)]dr{y). (21) 



Here a £W^, D is a real and positive semidefinite d x d matrix, z > and r is a probability 
measure on \ {0} s.t. its Fourier transform is entire analytic. The well-known interpretation 
of a is the drift vector, D determines the diffusion part, whereas z and r give the frequency and 
distribution of jumps, respectively. 

ip is called the generator of the (distribution of the) Levy process and the function 
k — > ijj{ik) is called the symbol of ip. 

4 graph expansion of the generalized OU process 



This section will be reserved to the graph applications to equations of type ([T]). We recall a 
graphical representation of the generalized OU process by a special type of graphs called rooted 
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trees. Let us start with some useful tools which help the reader to understand our formalism. 
Let G{t, x) be the Green function which satisfies : 



(22) 

G{t, x) = 0, t < ^ ^ 

Here 5{x) is the Dirac distribution. 

Proposition 4.1. The analytic solution of the generalized OU equation ([7P is given by: 

Xt = G* r/(x) + G * fix) (23) 
Proof. From the generalized OU equation we obtain 

(1+ m)Xi(x)= r/(x) (24) 
Now the result hold by applying the Green function given by equation 



Definition 4.2. A tree is a graph in which any two vertices are connected by exactly one simple 
path. In other words, any connected graph without cycles is a tree^ A tree is called a rooted tree 
if one vertex has been designated the root, in which case the edges have a natural orientation, 
towards or away from the root. 

Each vertex of a tree T connected with only one edge is called leaf. 

In our model we consider a rooted tree T with two types of leaves. The set of all trees with 
two types of leaves and i inner vertices will be denoted by T(i). 

Now it's of a great importance to ask about a method which gives an analytic value to each 
rooted tree T G T{i). The following algorithm named as "Feynman rule" gives a numerical 
value to each rooted tree with two types of leaves: 

Definition 4.3. Let TZ{T, rf) be the analytic value of a rooted tree T G T(i), then TZ is obtained 
as follows: 

- Assign r = [t, x) to the root of the tree T. 

- For each edge of the tree T multiply by G, where G is the green function given by equation 



For each leaf of type one multiply by the noise rj. 

For each leaf of type two multiply by the initial condition f . 

Integrate with respect to the Lebesgue measure dx. 



Definitions (14. 2p and (|4.3p can be summarized by the following table and graph: 

Theorem 4.4. Let Tj S 7~(i), j G N. The solution of the generalized OU equation ([7]) is given 
by a sum over all rooted trees that are evaluated according to the rule given in definition ^4-3^ > 
i.e. 

Xt= '^(T„r^) (25) 

TjeTi'i) 



^The various kinds of trees used as data structures in computer science are not really trees in this 
sense, but rather, types of ordered directed trees. 
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root 


inner vertices 


leaf 1 


leaf 2 


edge 


X 


• 


o 







Table 1: Different types of vertices. 



X- 




Figure 1: Construction of a rooted tree with two types of leaves and one inner vertex. 



Proof. From proposition ()4.ip one thus get: 

Xt = G* r]{x) + G* f{x) 
= X-® +X— 

= 7^(Tl, ?7) + 7^(^2, r?) (26) 



The graphical representations introduced previously can be applied for more complicated OU 
equations, e.g by replacing the right hand side of equation ([T]) by — mX^ + -q, in this case our 
graph formalism still true and the rooted trees will contain more inner vertices, hence in the 
rule given by definition (14. 3p we will integrate over all vertices. 
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